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Abstract 

White dwarfs and neutron stars are stellar objects with masses comparable to that of our 
sun. However, as the endpoint stages of stellar evolution, these objects do not sustain any 
thermonuclear burning and therefore can no longer support the gravitational load of their own 
mass by generating thermal pressure. Rather, matter in their interiors is compressed to much 
higher densities than commonly found in normal stars, and pressure is created by degenerate 
fermion kinetic energy and particle interactions. As a result, white dwarfs and neutron stars 
offer unique cosmic laboratories for studying matter at very high densities. In this review we 
discuss the basic properties of condensed matter at extreme densities and summarize the extent 
to which these properties can be examined by observations of compact objects. 



I. Introduction 



Astronomical phenomena provide many examples where matter exists in extreme conditions not 
found in terrestrial environments. One example is the high density of degenerate matter in "compact 
objects" - the relics of stars that have ceased burning thermonuclear fuel, and thereby no longer 
generate thermal pressure to support themselves against gravitational collapse. By contracting 
appreciably from their original sizes, the interiors of compact objects reach sufficiently high densities 
to produce nonthermal pressure via degenerate fermion pressure and particle interactions. Compact 
objects provide cosmic laboratories for studying the properties of matter at high densities. 



Firm observational evidence and well-founded theoretical understanding both exist for two 
classes of compact objects which support themselves against collapse by cold, degenerate fermion 
pressure: white dwarfs, whose interiors resemble a very dense solid, with an ion lattice sur- 
rounded by degenerate electrons, and neutron stars, whose cores resemble a giant atomic nucleus 
- a mixture of interacting nucleons and electrons, and possibly other elementary particles and con- 
densates. White dwarfs are supported by the pressure of degenerate electrons, while neutron stars 
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are supported by pressure due to a combination of nucleon degeneracy and nuclear interactions. 
These unique states of matter are achieved by significant compression of stellar material. Table || 
compares the principal physical quantities of a typical white dwarf and neutron star with those of 
the sunQ. 



Table 1: Parameters for the sun and a typical white dwarf and neutron star. 



object 


mass 


radius 


mean density 


mean pressure 


GM/Rc^ 




(Mq 


(km) 


(gm cm^^ ) 


(dyne cm~^ ) 




Sun 


1 


~ 7 X 10^ 


~ 1 


-IQl^ 


- 10-^ 


White Dwarf 


< 1.4 


~ 5 X 10^ 


~ 10^ 


~ 1024 


~ 10-4 


Neutron Star 


1-3 


~ 10 


-10^4 


~ 1034 


~ 10-1 



° - Mq = 1.989 X 10 gm = one solar mass. 

^ - This ratio measures the importance of relativistic gravitation, i.e., general relativity. 



Condensed matter in compact objects spans an enormous range of densities, which we loosely 
refer to as "high densities". These extend from about 7 gm cm-^ (e.g., the density of terrestrial 
IgFe), at the surface of a cold neutron star or white dwarf, to as much as p ~ 10^^ gm cm"^ , 
several times the density in atomic nuclei, in the cores of neutron stars. Matter at the various 
densities found in compact objects exhibits a variety of novel properties. Electromagnetic, strong, 
and weak interactions all play an important role in determining the character of compact objects. 
Since these objects are bound by gravity, they are a meeting point of all four of the fundamental 
forces of nature. Correspondingly, the astrophysics of white dwarfs and neutron stars incorporates 
a wide variety of physics including nuclear, particle, solid state and gravitation physics, to name a 
few areas. 

In this review we briefly survey the theory of condensed matter at high densities in compact 
objects and illustrate how the basic theory is tested through astronomical observations. Since we 
must cover fifteen orders of magnitude in density, our presentation is at most introductory in nature, 
and we encourage the interested reader to pursue the cited references. A detailed introduction to 
the physics of high density matter and compact objects can be found in the textbook Black Holes, 
White Dwarfs and Neutron Stars: the Physics of Compact Objects, by Shapiro and Teukolsky |]l|. 

We begin by considering the fundamental nature of cold (T = 0) high density matter in Sec- 



tion |L|, and distinguish between different regimes of high density. In Section |III.| we connect these 
microscopic properties with the fundamental macroscopic parameters of a compact object through 
the hydrostatic equilibrium dependence of mass and radius on central density. We summarize the 



fundamental predictions regarding the structure of white dwarfs and neutron stars. In Section IV. 



V. we examine how observations of white dwarfs and neutron stars can be used to probe the 



properties high density matter. We briefly discuss the perturbative effects of a finite temperature 



in Section VI 



^Throughout this review we will be using units which are the standard in astrophysical research: cgs for microscopic 
properties of matter and solar units (denoted by ©) for macroscopic properties of astronomical objects. 
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II. The Cold, High Density Equation of State 



The pressure and energies in white dwarfs and neutron stars are nonthermal; thermal effects due 
to a finite-temperature can be treated as a perturbation. We may therefore treat high density 
matter as having a zero temperature to very good approximation. The equation of state (EoS) 
of the matter then reduces to a single-parameter function, P{po) and p{po), or P{p), where P is 
the pressure, po is the rest-mass density and p is the total mass-energy density which accounts for 
internal (possibly relativistic) particle energies as well as rest-mass energy. 

There exist two main regimes of high density, distinguished as follows. As long as all nucleons 
are confined to nuclei, their contribution to the total pressure is negligible compared to that of 
the degenerate electrons. At some threshold density, Pn-drip, it becomes favorable for the nuclei to 
disintegrate, i.e., neutrons "drip" out of the nuclei and form a "nucleon gas". The standard EoS of 
Baym, Pethick and Sutherland [|| suggests that Pn-drip ~ 4 x 10^^ gm cm~^ . We may distinguish 
between the EoS for p < pn-drip which characterizes matter in white dwarfs and in the outermost 
layers of neutron stars, and p > Pn-drip^ which describes matter in the interior of neutron stars. 

A. The Equation of State below Neutron Drip Density: p ^ 4 x 10^^ gm cm~^ 

In matter below the neutron drip density the ions provide a Coulomb lattice of point-like charges, 
which is (to good approximation) independent of the properties of the surrounding electrons. The 
EoS of such matter is governed mainly by the electron gas. To lowest approximation, we may treat 
the electrons as an ideal fermion gas, incorporating some Coulomb corrections at relatively low 
density {p < 10^ gm cm~^ ) and corrections due to inverse /3— decay just below the neutron drip 
density (10^ gm cm~^ < P < Pn-drip)- The EoS of condensed matter below neutron drip density 
is well understood. The standard equations for cold, degenerate matter in white dwarfs (helium, 
carbon, oxygen and, possibly iron dominated models) have been derived by Chandrasekar and 
Salpeter Q, whereas models for equilibrated matter^, are based on the works of Dirac P and 
Feynman, Metropolis and Teller Q for p < 10^ gm cm~^ , and e.g., of Harrison and Wheeler [Q] 
and Baym, Pethick and Sutherland Q for higher densities. 

1. The Ideal Fermion Gas 

For almost the entire range of high densities the electrostatic energy associated with the structure 
of matter is much smaller than the Fermi energies. Consequently, Coulomb forces are generally 
negligible in a first order treatment of the high density EoS. The electron component of high density 
matter can therefore be described by a cold, single species gas of noninteracting fermions. At zero 
temperature the fermions fill all the states with momentum p < pf and none of the states with 
p > Pf, where pi? is the Fermi momentum. The corresponding Fermi energy of the particle species 

^The equilibrium isotope of matter is the nucleus of highest binding energy per nucleon. At low densities this 
isotope is normal 26^6, but as density increases so do the atomic mass and neutron to proton ratio 
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IS 

Ep ^ [{ppcf + {mc'fy^' , (1) 
where c = 3 x 10^'^ cm sec^^ is the speed of hght in vacuum and m is the fermion rest mass. 

For electrons, their number density, ng, is directly related to their Fermi momentum, pp^e by 
integrating over all occupied phase space {h = 6.63 x W~^^ ergs sec is Plank's constant): 

He = / ne{p)(i) = 2-. / 47r/45 = • (2) 



/o Jo 3/i3 

The factor of 2 in the second equation arises from the electron spin degeneracy. 

The pressure the electrons supply is calculated through the mean momentum flux of the electron 

gas, ^ ^ ^ 

where Xg = pF,e/'^eC is the electron "relativity parameter", Ag = h/{2Trmec) is the electron Comp- 
ton wavelength, Ve = PeC^/Ee is the electron velocity and 

^(X) = {x(l + x2)l/2(2^2/3 _ 1) + In + (1 + ^2)1/2] I _ 

The mass-energy density of the free electrons is also uniquely related to the Fermi momentum 

as 

o2^2 1 i' 2 



Ee{p)n^{p)dp = -3 / (p^c^ + (mc")2)^/^47rp> = ^xC^^e) ; (5) 
where 

= ^ + x2)i/2(i + 2x2) - In [x + (1 + x^ fl'^'\ } . (6) 

However, even when the degenerate electrons contribute most of the pressure, the mass-energy 
density is dominated by the rest mass of the ions, which are very nonrelativistic at these densities. 
Thus, the density of the matter may be simply taken as 

nerriB 

P = Po = -^rr- ' (7) 

J- e 

where Ye is the mean number of electrons per nucleon and is the mean nucleon mass. In the 
case of white dwarfs we may set Yg = Z/A = 0.5 {Z =atomic number, A =atomic weight) which is 
appropriate for fully ionized helium, carbon or oxygen, the most abundant constituents in a white 
dwarf, and uib = 1-66 x 10"^'* gm. Combining Eqs. (|2|), (|3[) and (^) provides the basic EoS of 
electron-pressure dominated high-density condensed matter in this regime. 

There exist opposite limits to the cold fermion gas equation of state: the low density, nonrel- 
ativistic (x <^ 1) and the high density, extremely relativistic (x ^ 1) limits. From Eqs. (^) and 
(^) we find that in cold matter with Y^ = 0.5 the electron relativity parameter satisfies Xe = 1 at 
p « 10^ gm cm~^ , which therefore marks the transition density between a nonrelativistic (NR) 
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and extremely relativistic (ER) electron gas. It is convenient to write the EoS in the two hmits in 
a polytropic form, 

P{p) = Kp^ (8) 

where 

NR < 1, p < lO^^gm cm-3 ; y = -, K = 1.0036 x 10^^^^^/^ ; (9) 

ER Xe » 1, p » lO^gm cm-3 : T = K = 1.2435 x IQ^^Y^/^ . (10) 

3 

The constant K is in cgs units, yielding a pressure in dyne cm^^ for a density in gm cm^'^ . Note 
that in this approximation the composition of the matter enters only through Yf>. Correspondingly, 
helium, carbon and oxygen, which all have Yg = 0.5 have identical ideal equations of state, slightly 
stiffer than that of iron {Y(, = 26/56 ~ 0.43). Fully equilibrated matter (often referred to as 
"catalyzed") has a Yf, that decreases with density and is therefore softer than matter composed of 
a single element. 

The free, degenerate electron pressure EoS outlined here is a good approximation for the equa- 
tion of state below neutron-drip density. It was employed by Chandrasekar in his pioneering 
analysis of equilibrium white dwarfs Q, for which he received the Nobel prize in 1983. More ex- 
act treatments were introduced in later years, which included the two main required corrections - 
electrostatic effects at low densities, and neutronization (or inverse /3— decay) at higher densities. 



2. Electrostatic Corrections to the Cold Equation of State: p <^ 10^ gm cm ^ 

There exists a net electrostatic correction to the ideal equation of state due to the fact that the 
local distribution of charge is very nonuniform. The fact that positive charge is concentrated in 
point-like ions causes the average electron-ion separation to be smaller than the average distance 
between electrons. The net electrostatic potential felt by the electrons is thus an attractive one, 
which effectively reduces the pressure for a given density. 

Electrostatic corrections to the cold equation state are mostly important at relatively low den- 
sities. Electrostatic energies are inversely proportional to the average separation between particles, 

— 1/3 

< r > which is naturally proportional to rig . The relative importance of electrostatic energy, Ec, 
between a degenerate, nonrelativistic electron and an ion of charge can be estimated through 

|£ . . „-/3 (u) 

where E'p^ = oc n^^'^ is the Fermi kinetic energy of the nonrelativistic electrons. Unlike the 
case of hot matter (where the mean electron kinetic energy is ~ ksT), the relative importance of 

— 1/3 

electrostatic corrections decreases with density as ne 

A rough estimate of the quantitative electrostatic correction can be performed by using the 
Wigner-Seitz approximation, which describes the lattice as neutral sphere with a central point- 
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like ion and an ambient uniform electron ga^. For a carbon or oxygen lattice, the correction to 
the pressure is typically of the order of a few percent. At lower densities the electron distribu- 
tion deviates from uniformity, and more sophisticated approaches, namely the Thomas- Fermi and 
Thomas- Fermi- Dime models, must be invoked [Q, 

3. Neutronization Corrections to the Cold Equation of State: 10^ gm cm p ^ 4 x 

10^1 gm cm^^ 

Stable high density matter must be in chemical equilibrium to all types of reactions, including the 
weak interactions which drive P decay and electron capture ("inverse /?— decay"): 

n^p + e + Pe , p + e^n + i/g, (12) 

where n and p denote a neutron and a proton, respectively, e denotes and electron and fe(^'e) 
denote an electron neutrino (anti- neutrino). If the matter's composition is out of /5— equilibrium, 
it will adjust through /3— decays or electron capture. Both types of reactions change the electron 
per nucleon fraction, Yg, and thus affect the EoS (Eqs. (p|-px|)). 

In cold white dwarfs and neutron stars the weakly interacting neutrinos freely escape the system: 
a zero neutrino abundance implies a zero neutrino chemical potential. The condition of chemical 
equilibrium is then stated as: 

fJ-n = fJ-p + fJ-e , (13) 

where fj,x denotes the chemical potential of species x. The condition of chemical equilibrium is the 
fundamental origin of the stable existence of neutrons in nuclei and in uniform n — p — e matter. 
For free, single, particles the chemical potential is identical to the rest mass. The masses in MeV 
of the three elementary particles of Eq. ( p!3[ ) are m„ = 939.6, rup = 938.3 and mg = 0.511 (1 MeV 
is equivalent to ~ 1.78 x 10"'^^ gni). It is energetically allowed to have a free neutron decay in the 
reaction n^p+e + Pe {^n — f^p — ^nrie ~ 0.8 MeV). Indeed, the life-time of a free neutron is only 
about 1000 sec before it undergoes a /3-decay. By contrast, in a cold, degenerate, noninteracting 
n — p — e gas neutron decay can be "blocked" : since the protons and electrons must obey the Pauli 
principle, the decays will be suppressed if the energy available to the newly formed electron and 
proton is insufficient to place them above their respective Fermi levels. At such high densities, the 
equilibrium state of the gas includes a finite fraction of neutrons. 

In bulk matter, the electron chemical potential, which is equal to the electron Fermi energy, 
rises with electron number density. Hence, maintaining chemical equilibrium (Eq. |l^)) may require 
some protons to capture electrons and covert to neutrons. In matter below the neutron drip density 
these conversions occur in nuclei, so their neutron fraction increases - they "neutronize" . The result 
is a net decrease in the electron abundance (Ye) at high densities, which lowers the pressure and 
"softens" the EoS. 

The fact that "neutronization" becomes energetically favorable at densities below the neutron 
drip must be taken into account when formulating an exact EoS for this range. In matter composed 

^We note that such an approximation is actually better suited for white dwarfs than for laboratory solids, where 
the electron distribution is much more nonuniform. 
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of nuclei and free electrons, /?— decay is limited to relatively high densities (p > 10^ gm cm~^ ), and 
the exact threshold for the onset of neutronization depends on the nature of the nucleus, since the 
nuclear interactions must be taken into account in addition to Fermi kinetic energies. For example, 
the reaction g^C + e — >P B requires the nucleus to absorb an energy of ~ 14 MeV, which an 
electron can supply if the density of a pure carbon lattice reaches ~ 3 x 10^'' gm cm~^ . The EoS 
of matter close to the neutron drip density is thus dependent on its chemical composition even for 
elements of equal Y^, such as helium and carbon |^. 

B. The Equation of State Above Neutron Drip Density: p 4 x 10^^ gm cm~^ 

As neutronization proceeds, the nuclei become increasingly neutron-rich. The so-called "tensor- 
force" component of the nuclear interactions causes like nucleons to repel one another, so that 
the binding energy of a neutron rich nucleus is smaller than in one where Z/A = 0.5. Fully 
equilibrated matter, reaches the last stable isotope ^^'^Kr (Y^ 0.31) just before the density of 
p = Pn-drip ~ 4.3 X 10^^ gm cm~^ . At higher densities it becomes favorable for some neutrons to 
"drip" out of their parent nuclei. The nucleon component of the matter can no longer be confined 
to point-like objects and the nuclei begin to dissolve. As density is increased, a larger fraction of 
the nucleons exists as "free" particles, outside the nuclei. This is a gradual transition, until the 
matter approaches the density of atomic nuclei, where all nuclei have essentially dissolved, and 
the distribution of nucleons becomes uniform. The EoS above the neutron drip density must be 
formulated by consistently including the effects of nuclear physics, which become the governing 
component of the properties of matter as the density of atomic nuclei, p^^c ~ 2.8 x 10^^ gm/cm^ is 
approached. 

C. Subnuclear Densities: 4 x 10^^ gm cm~^ ^ P ^ 2.8 x 10^^ gm cm~^ 

An analysis of cold matter in the range Pn-drip ^ P ^ Pnuc is complicated by requiring chemical 
equilibrium between nucleons inside the nuclei and those that have dripped outside. One must 
account for the effects of the surrounding gas of free nucleons on the nuclei, as well as other effects 
such nuclei surface and Coulomb energies. In addition, the nuclei are expected to be very neutron 
rich, deviating from the Z/A = 0.4 — 0.5 found in terrestrial nuclei. Nonetheless, the properties 
of matter in this range of densities can still be derived by a natural extrapolation from ordinary 
nuclei, and, indeed, the EoS for such matter is believed to be well understood. The principal 
studies of matter (e.g., Ref. [^) are based on the nuclear liquid-drop model, and are suitable for 
most applications regarding neutron star structure. For problems where a more accurate description 
of the neutron star inner crust is required, attention must be given to lattice effects, as the nuclei 
and free nucleons arrange in a distinct spatial structure (where the nuclei settle into bubbles, slabs 
or rods, depending on density [^). It is expected that at p « ^Pnuc all nuclei will have dissolved so 
that the matter is completely uniform. 
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D. Supernuclear Densities: p ^ 2.8 x 10^'' gm cm ^ 

The upper end of the high density regime is referred, to cis "supernucles-r" where p ^ Pnuc- 

As 

matter is compressed to such densities the EoS becomes gradually dominated by the degenerate 
nucleons and the nucleon interactions. 

It is instructive to begin by considering an ideal uniform mixture of neutrons, protons and 
electrons. For any given nucleon number density, titv, solving the equilibrium composition n„,np 
and Ue (the neutron, proton and electron number densities, respectively) requires three equations. 
The first is baryon number conservation, 

UN = Un + Up , (14) 

and the other two arise from imposing chemical equilibrium (Eq. |l3|) and charge neutrality. The 
condition for charge neutrality is 

rip = rie . (15) 

For a noninteracting mixture of nucleons and electrons all chemical potentials are simply the 
Fermi energies (Eq.^)). Electrons are extremely relativistic at nuclear densities so their electron 
chemical potential is 

lie - PF,eC = hc{?,^^Y,pBf'^ « 100 (^) MeV , (16) 

V0.03/ \nnucJ 

where Unuc = 0.16 fm~^ is the number density of nuclei at the saturation density (and 1 fm = 
10^^'^ cm). Even if the electron fraction per baryon is only a few percent, the electron chemical 
potential exceeds the mass difference between neutrons and protons by two orders of magnitude. 
The only way to maintain a finite electron fraction (which is required to balance the protons for 
charge neutrality) and satisfy chemical equilibrium is by having a significantly larger neutron than 
proton fraction (note that for nonrelativistic fermions Ep oc n^/^ and for extremely relativistic 
fermions Ep oc n^^^). Equilibrium matter at nuclear densities must be very neutron dominated, 
and objects composed of such matter are thus "neutron stars". For an noninteracting gas the ratio 
of neutrons to protons in equilibrium must be about 8:1, which is also representative of more 
realistic models of supernuclear densities. 

The noninteracting gas approximation is not reliable for deriving the EoS at supernuclear 
densities. Unlike electrostatic perturbations, nucleon-nucleon interactions are not negligible, and 
the interaction energies are comparable to the Fermi energies of the degenerate nucleons (electrons 
do not feel the strong interaction, and may still be treated as noninteracting). Modeling of the 
nucleon-nucleon interaction is one of the longest-standing problems in nuclear physics, still only 
partially solved. Profound difficulties exist due to the absence of a comprehensive theory of the 
interactions and the difficulty of obtaining experimental data for p > Pnuc- Further complications 
arise due to the fact the Fermi and interaction energies at p ~ 2 — 3pnuc reach a sizable fraction 
of the rest mass, and relativistic effects must be taken into account as well. It is not yet possible 
to apply quantum chromodynamics (QCD), the fundamental theory of strong interactions, to the 
many-body nuclear domain at p ~ Pnuc- Instead, the most useful approaches are still based on 
phenomenological potential formalisms, and many-body Schrodinger-like systems of equations Q. 
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The EoS of supernuclear matter remains to date a field of active research. Current approaches 



include variational methods based on deduced two and three nucleon interactions ( , see ||lT| for 
a recent study), and relativistic mean field approximations (|jl2|, see [^] for a review). The models 
are based on fitting parameters to reproduce the empirically determined properties of finite nuclei 
with n = Unuc- A rough approximation of the properties of the nucleon component of supernuclear 



matter can be obtained through the effective, nonrelativistic model of Ref. |14|, 

p{n, x) = Egas{n, x) + an^ + - Upf + cn^^^ . (17) 

The first term Egas is the mass-energy density of a noninteracting gas of nucleons at density n and 
composition nn,np. The coefficient a is negative, representing the long range attractive component 
of the inter-nucleon force, while c is positive, representing the short-range repulsive component. 
The power 5 is larger than unity, so that short-range repulsion dominates at high density. The 
symmetry term includes the positive coefficient b which describes the "tensor-force" that repels like 
nucleons, and its contribution is therefore minimized at for symmetric nuclear matter (n„ = Up). 



The values of a, b, c and 6 are derived by requiring Eq. (17) to reproduce the assumed properties of 
symmetric matter at hn = nnuc- See Refs. |l5| for typical (model dependent) values of these 
parameters. 

Finally, we note that it is quite possible that other particles, besides neutrons, protons and 
electrons, coexist in stable equilibrium at supernuclear matter. The most obvious example is the 
muon, which is a lepton similar to the electron, but has a rest mass of ~ 105 MeV. The 
equilibrium condition for the muons will simply be /x^ = fie, so from Eq. ([l^ ) it is evident that 
at densities n ^ n-nuc it is energetically favorable to convert some electrons into muons through 
the weak interactions. At densities of n '^nnuc it possible that other exotic particles will appear, 
such as hyperons (baryons that are heavier than nucleons), Bose-Einstein condensates of mesons 
(i.e., pions or kaons), or even conversion of the nucleons into an uniform mixture of quarks (see 
Refs. |l3| , H^, ^] for recent reviews of the possible presence and roles of such particles in neutron 
stars). 



E. Basic Properties of High Density Matter 

In Table |2| we present an EoS of cold, fully catalyzed matter ranging from normal iron at zero 
pressure to supernuclear densities. The EoS in the subnuclear regime is compiled from Refs. [§, 
|6|, ^ and are "standard" for studying catalyzed high density matter. For realistic white dwarf 
models, which are presumably composed mostly of helium or carbon and oxygen, and never reached 
sufficiently high temperatures to catalyze their nuclei, a slightly different EoS must be used, based 
on a single species; see Ref. 0]. We also tabulate one "state-of-the-art" EoS for the supernuclear 



range |11| 



The EoS P{p) is plotted in Figure |T|a. It is especially instructive to examine corresponding 
values of the adiabatic index, F = dlnP/dlnp, plotted in Figure |l|b, along with the sound speed 
(cs = y/dPjdp). Note that around p ~ 10^ gm cm~^ the EoS transforms, as expected, from a 
nonrelativistic F = 5/3 polytrope to a relativistic F = 4/3 one as the Fermi energy of the electrons, 
which dominate the pressure, gradually becomes relativistic. There is a sharp drop in the adiabatic 
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Table 2: 


The equat 


ion of state 


° of cold, 


catalyzed high density matter 


P 


Ub 


P 


r 


y 


iv 


Ca 


gm cm~^ 




dyne cm~^ 




dyne cm 


ivie V 


cm sec ^ 


l.OOE+01 


6.02E+24 


5.13E+11 


7.03E+00 






K CTTP 1 nc^ 


2.00E+01 


1.21E+25 


2.59E+13 


3.09E+00 


7.o2i^+lo 






5.00E+01 


3.01E+25 


2.21E+14 


2.50E+00 
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2.16E+22 


1.59E+00 
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1.20E+30 
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3.01E+30 


2.63E+23 
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" - shortened notation: 2.00E + 10 is 2.00 X IQl". 

- The EoS up to p ~ pnuc = 2.8 X lOli is based on Table 5 of Ref. |^ [© 1971. The American Astronomical Society]. The 
supernuclear EoS is adapted from Ref. [hl| . 

- the APR | |l]| EoS is not relativistic, and at extremely high densities has an unphysical superluminal (greater than light 
speed) speed of sound. 
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Figure 1: The equation of state of cold, catalyzed high density matter, based on Refs. [P, ^, 
(a) the pressure density relation, (b) adiabatic index and sound speed. 



index around the neutron drip density, since at first the dripped neutrons contribute to mass-density 
but not to pressure. Only as matter approaches nuclear densities does the nucleon pressure become 
important, pushing the adiabatic index back up to values in the range 2 — 3. 

Along with the EoS, we also examine some thermodynamic quantities of cold high density mat- 
ter. In Table m we also list the adiabatic index, sound speed, bulk modulus Y, and incompressibility, 
K, defined as 

Y = n^ and K = 9^ = 9n^ . (18) 
dn dn dn'^ 

Note that the bulk modulus (which has units of pressure) is the reciprocal of the quantity usually 
defined as the compressibility {x = ^^^)i while the "incompressibility" (units of energy) is more 
commonly used in nuclear physics applications (and is indeed measurable in nuclei) . In Table ^ we 
list typical values of T, c^, Y, K and the specific heat capacity, c^, for condensed matter found in 
white dwarfs and neutron stars. For the purpose of comparison, we also give values for the sun, 
where the pressure roughly follows an ideal Maxwell-Boltzmann law P ~ nksT. Note how the 
extreme conditions of high densities leads to properties that are very different than those found for 
terrestrial materials. For example, the speed of sound in a neutron star is expected to be several 
tens of percent of the speed of light! 
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Table 3: Typical values of thermodynamic quantities in the sun, white dwarfs and neutron stars . 

Object r Y K c„ 

(cm sec^^ ) (dyne cm~^ ) (MeV) ergs deg~^ gm"""^ 
~~S^ 5/3 ^^^^ ~ 10"^ ~ 10"'' ~ f (l//i + l/mp)fcs 

White Dwarf (C/O) 4/3 ~ 5 x 10** ~ 10^* - 10 ffcs//^ 

Neutron Star 2-3 ~ lO"' ~ 10^^ ~ 200 ~ lO'^fca/mp 1^ 

t - for an ideal nonrelativistic neutron gas and a central temperature of ~ 10^ °K — > ksT/Ep ~ 10"'^. 



III. White Dwarf and Neutron Star Structure 



Compact objects are self-gravitating equilibria with a mass comparable to that of the sun, 
1 Mq = 1.989 X 10^^ grams. Both white dwarfs and neutron stars are centrally condensed ob- 
jects - most of their mass is located in a high density core, which is limited to a fraction of the 
volume; furthermore, the radius of the objects decreases with increasing mass. These traits are 
characteristic of a configuration supported by degenerate-fermion pressure. 



A. Construction of a Self-gravitating, Equilibrium Star 

The equilibrium structure of a self- gravitating object is derived from the equations of hydrostatic 
equilibrium. The simplest case is that of a spherical, nonrotating, static configuration, where 
for a given EoS all macroscopic properties are parameterized by a single parameter, for example, 
the central density. In the case of compact objects, the gravitational fields are strong enough 
that calculations must be performed in the context of general relativistic (rather than Newtonian) 
gravity. The fundamental equation of hydrostatic equilibrium in its general relativistic form has 
been derived by Tolman [^] and Oppenheimer and Volkoff [^, and is known as the "TOV" 
equation: 

dPjr) ^ Gm{r)p{r) / ^ P{r) \ / ^ A7:r^P{r) \ / 2Gm(r) \-^ 
dr \ c^p{r) J \ c^m{r) I \ c^r J 

This equation simply states that at any radial distance r, the gravitational pull by the mass interior 
to r, m(r), is balanced by the gradient of the pressure P{r); G = 6.67 x 10~® cm^ gm~^sec~^ is 
the gravitational constant. Note that the first term on the right hand side is the only term in the 
nonrelativistic case, where P/{p(P') <C 1 and 2Gm/ {re?) <C 1. The second and third factors arise 
from the pressure being a form of energy density ("regeneration of pressure" effect), while the last 
term includes the correction due to the curvature of space in the strong gravitational field of the 
star. 
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B. Stable Configurations of High Density Matter 



Given an EoS, Eq. (19) can be integrated simultaneously with the mass equation, 



dm{r) 
dr 



47rr^p(r) 



(20) 



to determine the entire profile of the object. The actual integration is typically solved by numerical 
means, assuming a central density p{r = 0) = pc, and integrating outward from the center until 
reaching the surface where P(r) = 0, which identifies r = Ras the radius of the star and M{R) as its 
mass. We emphasize that the TOV equations provide the gravitational mass (or total mass-energy) 
of the star, which includes the effect of the gravitational binding energy, Egb ~ GM'^/R, as well 
as the internal and rest mass energy of the stellar constituents. Note that in the case of neutron 
stars, {GM'^/R)/{Mc^) ~ 10-20%; indeed, the gravitational mass of a neutron star is measurably 
lower than the total rest mass of its constituents. 



2.5 



1.5 




^■^ // cold-catalyzed matter 



10' 10^ 10"' ^lo" 
p^, (gm/cm ) 



lo" 10" 



Figure 2: Mass vs. central density of a cold, carbon star (dashed line) and of a star composed 
of cold catalyzed matter (solid line). 



Varying pc produces a sequence of models for the given EoS, yielding M{pc) and R{pc) along the 
sequence. The resulting sequences for stars composed of pure carbon (based on |]l8|) and of cold, 
catalyzed condensed matter (with the EoS tabulated in table 0) are both shown in Figure I (note 
that carbon dominated matter beyond the neutron drip density is unphysical and therefore omitted). 
The most distinct feature about the equilibrium sequences is the existence of local maxima in the 
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White Dwarf 



Neutron Star 




Figure 3: Schematic illustration of the structure of a typical white dwarf (left) and a neutron star 
(right). 



M{pc) curves. A hydrostatic equilibrium configuration is dynamically unstable to catastrophic 
gravitational collapse if dM/dpc < 0, since a radial perturbation would cause it to collapse on itself 
Q. Therefore stable carbon white dwarfs exist only if the central density is < 6 x 10^ gm cm~^ , 
while cold, catalyzed, matter has two distinct regimes of stable configurations: one with central 
densities below ~ 10^'^ gm cm~^ (equivalent to white dwarfs), and one where the central density lies 
roughly in the range 10^^ < Pc < 10^^ gm cm~^ (neutron stars). Configurations where pc ^ Pn-drip 
are unstable, and cannot be found in nature. This situation is a direct consequence of the nature 
of the EoS as shown in Figure |l|. Quantitatively, it can be shown that a hydrostatic configuration 
is stable only if its mass averaged adiabatic index, F, satisfies Q 

- 4 2GM , , 

It is evident therefore, that the effect of neutronization and neutron drip, which cause F to drop 
below 4/3, is to separate astrophysical objects with high density matter into two distinct classes of 
stable configurations, which are shown schematically in Figure |^. 



1. White Dwarfs 



White dwarfs are composed of matter below the neutronization density. The astrophysical scenario 
which creates white dwarfs, i.e., quasi-static contraction of a progenitor star at the end point 
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of thermonuclear burning, does not allow matter to reach high enough temperatures to achieve 
equilibrium composition through further burning. Accordingly, the composition of white dwarfs is 
mostly dominated by its nuclear ashes. In the standard scenario white dwarfs are composed mostly 
of carbon and oxygen, but there is observational evidence that helium white dwarfs (and perhaps 
even iron-core white dwarfs - see below) exist as well. The pressure throughout the star is due to 
degenerate electrons, which are nonrelativistic in the outer layers, but are relativistic in the interior 
of the more massive stars with M ~ Mq. 

By approximating the EoS as a Newtonian polytrope (Eq. (§)), Chandrasekar H] derived the 
basic features of white dwarfs: 



Pe^loVcm-3 r = |: R^l.l2xW'x( ^^\. ] ' {^] cm (22) 



1/6 ^Y^^5/6 
lO^gmcm^'^/ \0.5, 

1/2 / y N 5/2 



M ~ 0.496 X ( —rr-^ ^ 1 ( — 1 , 

V106gmcm-3y Vo.5/ 

/ o \ ^^^^ / Y \ 

Pc > lO^gm cm-3 r = I : i? ~ 3.35 x 10^ x ^ ^ — cm (23) 

^ 3 - VlOSgmcm-sy Vo.5/ ^ ' 

M ~ 1.46 X (^^^ Mq . 

White dwarfs are expected to have radii of the order of 10^ km - roughly the size of the Earth, or 
about one percent of that of the sun (which has a radius i?0 ~ 7 x 10^ km). 

As long as the EoS is approximated as a pure polytrope, so that electrostatic and neutronization 
corrections are ignored, helium, carbon and oxygen white dwarfs (Ye = 0.5) are all identical, while 
an iron-dominated white dwarf is only slightly different (1^ ~ 0.43). Cold, catalyzed matter, which 
is used for the mass sequence in Figure ^, has a significantly softer EoS (due to decreasing 
with increasing density), and its sequence lies, therefore, lower than those of stars composed of a 
single species. The most significant result in Eq. ( p3|) is that there is no dependence of M on the 
central density. The mass Mch — 1-46 Mq, known as the Chandrasekar mass, is the asymptotic 
value a white dwarf can reach if it achieves sufficiently high density, so that its entire structure is 
governed by relativistic fermions (in practice this mass cannot be reached, since the outermost layer 
have nonrelativistic electrons). Nonrotating white dwarfs cannot have a mass exceeding Mch^ and 
the precise limit on their mass is lower by several percent, due to neutronization at high densities 



(Section [n.B.3.) and general-relativistic effects. In spite of the negligible effect on the hydrostatic 



equilibrium profile of white dwarfs, general relativity is required for a full analysis of white dwarf 



stability, since an exact T = , star is unstable to gravitational collapse (Eq. |21 



2. Neutron Stars 



Neutron stars are composed mostly of matter at nuclear densities, including a core with supernuclear 
densities, and is topped by a thin crust at subnuclear densities. The crust, composed of cold, 
catalyzed matter, may also be divided into an inner part with pn-drip ^ P ^ Pnuc and an outer 
part where p < Pn-drip- 
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Immediately after the neutron was discovered (1932) Landau modeled a neutron star as a gas of 



noninteracting, degenerate neutrons 1 21 ] . He found that neutron stars would have a maximum mass 
of about IMq and a radius of several kilometers. Although this was a very crude approximation, 
it does suggest the correct orders of magnitude regarding the structure of these objects: a neutron 
star has a mass comparable to that of the sun compressed to the size of a medium city! If supported 
by degenerate fermion pressure, a neutron star cannot have a mass that exceeds the Chandrasekar 
limit by much, while its radius is inversely proportional to the mass of the pressure providing 
fermion. A neutron star should indeed have radius about m^/mg ~ 10^ times smaller that of a 
white dwarf. 
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Figure 4: The mass vs. radius relation for of a nonrotating neutron star for several supernuclear 
equations of state: Huber et al. (RHF8, [^2|), Pandharipande and Smith (MF, ||2^), Glendenning 
(GM24, [^], page 244), Friedman and Pandharipande (FPS, [p^]), Wiringa, Fiks, and Fabrocini 
(AU, |l^), and Akmal, Pandharipande and Ravenhall (APR, [|ll| ). 



A realistic treatment of neutron star structure requires general relativity [20|, whose effects are 
appreciable in this case. Furthermore, quantitative estimates must be based on a realistic EoS at 
all densities and especially a realistic model for the supernuclear regime where particle interactions 
are most important. In particular, these interactions are repulsive at short distances and oppose 
compression, thereby stiffening the EoS in comparison to a free particle gas. Indeed, plausible 
equations of state predict that Mmax{N S) ~ 2 Mq, (~ 2.2 Mq for Ref. |11]), whereas the star's 
radius will lie in the range i? ~ 10 — 15 km. As an example, we show in Figure ^ the M vs. R 
relation found for a few representative equations of state of supernuclear densities. The exact 
value of Mmax{N S) depends on the assumed EoS and therefore provides an integral measure of 
the properties of matter at supernuclear densities [|l7t 25 1. We also note that the existence of a 
maximum mass for neutron stars has important astrophysical implications, since it suggests that a 
larger mass cannot be sustained by cold pressure and must inevitably collapse to a black hole. 
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IV. Observations of Condensed Matter in White Dwarfs 



Over 2000 white dwarfs have been discovered to date, based on the spectroscopic properties of 



observed stars [^6|. While the thermal energy in white dwarfs is small compared to the Fermi 
energies of the electrons in the interiors, it is still sufficient to generate an observable surface 
luminosity for several billion years. The luminosity of a star is dependent on its radius, R, and its 
effective surface temperature, Tf,ff, according to 

L = AnaR^T^jf , (24) 

where a = 5.67 x 10~^ erg cm~^ sec~^ deg~^ is the Stephan-Boltzmann constant. The luminosities 
of white dwarfs are typically 10~^ to 10~^ of that of the sun, Lq 4 x 10^^ ergs sec~^, but their 
radii are also smaller than those of ordinary stars by a factor of ~ 100. As a result, white dwarfs 
have apparent surface temperature of several times 10^ °K , unusually high in comparison with 
most observable stellar objects, so they do appear very "white". The existence of white dwarfs was 
established spectroscopically (by determining the surface temperature of sources) as early as 1910. 

White dwarfs are believed to be remnants of stars with initial masses in the range 0.1 — 8 Mq, 
which are not massive enough to complete the thermonuclear burning process all the way to iron. 
Mostly they are mostly composed of carbon and oxygen, but in some cases, thermonuclear burning 
ceased before these elements were produced, and such stars are dominated be helium. 



A. Radii 



Observational determination of a white dwarf radius is straightforward if its flux, F, is measured 
and its distance from the Earth, D, is known: 

The stellar radius is derived only after the effective surface temperature is obtained spectroscopically 
p^ . Estimated radii of white dwarfs reside in the range of 0.007 — 0.013 i?0, which is consistent 
with an object supported by degenerate electron pressure (Section [[II.|B. ) and thereby confirms the 
basic nature of white dwarfs. 



B. Mass-Radius Relations 



The most significant test of the nature of matter in a white dwarf is obtained by comparing 
observed mass-radius relations with the theoretical predictions. Once the radius of a white dwarf 
is determined, details of the surface emission, namely effects of gravitational acceleration on line 
emission!^] (which scales as M/i?^), or gravitational redshift (which scales as M/i?), are then 
used to estimate the surface gravity. An independent (and usually more accurate) estimate of the 
mass can be obtained directly for white dwarfs in binary-star systems, through Kepler's third law 
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Table 4: Mass-radius relations of selected white dwarfs with la errors. [Adapted from Provencal 
et al. |27], © 1998. The American Astronomical Society] 



Object 



Mass 

(Mo) 



Radius 

(^0) 



Masses from binary motion 
SiriusB 1.000 ±0.016 0.0084 ± 0.0002 

Procyon B 0.604 ± 0.018 0.0096 ± 0.0004 
40EriB 0.501 ±0.011 0.0136 ± 0.0002 

Masses from surface acceleration 
EG 50 0.50 ±0.06 0.0104 ± 0.0006 

GD 140 0.79 ± 0.09 0.0085 ± 0.0005 

Masses from surface gravitational redshift 
CD-38 10980 0.74 ± 0.04 0.01245 ± 0.0004 
W485A 0.59 ± 0.04 0.0150 ± 0.001 



j33| |. Through these methods, masses of a few white dwarfs have been estimated with reasonable 
accuracy, and most seem to cluster around ~ 0.6 Mq [^1|. Some larger mass white dwarfs are also 
known, including the most famous Sirius B |34|, with M ~ 1 Mq. Note that these masses are 



smaller than Mch- it is believed that stellar evolution, and especially periods of mass loss, limit the 
masses of most white dwarfs to M < 1 Mq, even though the EoS could sustain somewhat higher 
masses. 

Some examples of observationally determined radii and mass for white dwarfs are presented in 
Table |[ Note that lar ger mass stars have smaller radii, as is expected (Eqs. (|2^p3|)). It is interesting 
to note that for some time, there was a nonnegligible discrepancy between the observations and 
theory, where observed radii seemed to be ~ 10 — 20% smaller than estimated by theory for carbon 



white dwarfs [29|. Only improved estimates of distances to several white dwarfs with the Hipparcos 



satellite, and better modeling of white dwarf atmospheres |35] have allowed this discrepancy to be 
mostly resolved [^] . A comparison of current observed mass-radii determinations with theoretical 
curves is presented in Figure ^, and in general, the fit is indeed very good: these new results 
seem to confirm that the composition of most white dwarfs is indeed dominated by medium weight 
elements (carbon and oxygen). However, they also imply that a small minority of white dwarfs 
do have relatively small radii, and may therefore contain iron cores, which presents an intriguing 
puzzle from the point of view of stellar evolution. 
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Figure 5: The mass-radius relation for white dwarfs. Sohd Hues labeled He, C, Mg, 



and Fe denote 
Models 



the zero-temperature mass-radius relation for a star composed of each element of Ref. [18| 
of white dwarfs with a hydrogen atmosphere [^] are for an effective atmosphere temperature of 
30,000°K (dotted hne) and 15,000°K and 8000 °K (dashed lines). Estimates for observed white 
dwarfs are shown with a la error. [From Provencal et al. [p7| ], (c) 1998. The American Astronomical 
Society] 




.006 



V. Observations of Condensed Matter in Neutron Stars 



Neutron stars are often identified observationally as very accurately pulsating sources, where the 
pulsation is attributed to the star's rotation. The period of pulsation limits the size of the source 
to no more than a few tens of kilometers, thus indicating the presence of a very compact object. 
Since the first discovery in 1967, close to 800 pulsars have been identified, over 700 of those in 
radio waves |36|. About thirty X-ray pulsars have also been observed, and a smaller number of 



nonpulsating X-ray sources of various types are also most likely to contain a neutron star 1 37 ] 



Theory suggests that neutron stars are formed when a massive star is disrupted in a "supernova" 
explosion. As the massive star evolves, matter in its core undergoes thermonuclear burning all the 
way to iron, which is the most tightly bound nucleus and therefore cannot burn further. When 
the iron core is massive enough it collapses under its own weight, until the collapse is halted when 
nuclear densities are reached. Most of the gravitational binding energy is released in the form of 
neutrinos, but some is transformed into an outgoing shock wave that expels the star's envelope and 
leaves behind the newly formed neutron star. Indeed, several dozen young neutron stars have been 
found in sites of past supernovae, confirming this scenario. 
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A. Neutron Star Masses 



The masses of over twenty neutron stars have been determined observationally through their grav- 
itational pull on a binary companion star. Within the errors of measurement, the masses of these 
stars appear to cluster in a narrow range of 1.35 it 0.1 Mq jS^. This is especially true of the 
masses of neutron stars in the four known binary-pulsar systems (where the companion object is 
also a neutron star), listed in Table |5|, where mass determinations are most accurate. Another 
example is the neutron star - white dwarf binary (once believed to be a double binary neutron 
star system) B2303+46. The properties of the orbit allow to limit the mass of the neutron star to 
1.2 < Mns < 1-4 Mq In particular, the Hulse- Taylor binary pulsar 1913+16 (the coordinates 
of its location in the sky) has been very accurately determined to have M1913+16 = 1.4411+0007 Mq. 
Clearly, a high density EoS must satisfy M^ax > 1-44 Mq to be consistent with this observation. 

Table 5: Mass determinations with errors of neutron stars in known binary pulsar systems. 
[Adapted from Thorsett and Chakrabarty @ 1999. The American Astronomical Society]. 





Median mass 


68% 


95% 


Star 


mass {Mq) 


central limits 


central limits 


J1518+4904 pulsar 


1.56 


+0.13/ - 0.44 




J1518+4904 companion 


1.05 


+0.45/ - 0.11 




B 1534+ 12 pulsar 


1.339 


+0.003 


+0.006 


B 1534+ 12 companion 


1.339 


+0.003 


+0.006 


B1913+16 pulsar 


1.4411 


+0.00035 


+0.0007 


B1913+16 companion 


1.3874 


+0.00035 


+0.0007 


B2127+11C pulsar 


1.349 


+0.040 


+0.080 


B2127+11C companion 


1.363 


+0.040 


+0.080 



This limit does not impose a serious constraint on realistic models for the supernuclear EoS and 
is satisfied by all but the softest equations. There has been much debate whether the narrow range 
of observed neutron star masses is evidence for the maximum mass being rather low {Mmax{N S) ~ 



1.5 Mq) [40|, or whether it is due to astrophysical effects, which could restrict the range of masses 
of observable stars. Recently, there is growing evidence that heavier neutron stars exist: one in the 
Vela X-1 binary is now estimated to have a mass of ~ 1.9 Mq [^1[, and the oscillations observed 
in some galactic X-ray sources seem to indicate that they include neutron stars with masses larger 
than 2 Mq |^]. If these estimates are confirmed with improved observations in the future, they 
will provide a much more stringent test that would be consistent only with the stiffer models for 
the supernuclear EoS. 



B. Radii 

The mass and radius relation of a nonrotating neutron star is uniquely defined for any given EoS 
through the solution of the TOV equations (Figure §). However, to date, observational methods for 
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estimating neutron star radii lack the accuracy required to critically differentiate between realistic 
equations of state. Hopefully, data from a new generation of X-ray satellites will provide substantial 
constraints on the compactness (M/R) of neutron stars (either through their thermal emission 
or from emission of material accreting upon them p^ ). 



C. Rotation periods 

The observed pulsation period of a neutron star is attributed to rotation. Most measured periods 
are in the range 0.25 — 1 sec, but a subclass of millisecond pulsars are also known to exist. The 
fastest spinning neutron star observed to date has a period of P = 1.56 msec (so that it spins faster 
than an egg-beater!). 

Rotation provides an additional centrifugal barrier which assists the internal pressure in sup- 
porting the gravitational load. The resulting configuration is then dependent on both central 
density and rotation period, and must be solved self-consistently including the effects of general 
relativity. The maximum mass a given equation of state can support increases with respect to the 
static (TOV) value as the star is allowed to rotate faster [^]. There must exist a lower limit on 
the rotation period (upper limit on the rotation rate), since too rapid rotation would cause the star 
to shed mass at the equator. The general trend is that a softer EoS predicts a smaller lower limit 
on the rotation period, since a more compact star is more tightly bound gravitationally and so is 
better suited to resist the centrifugal force. At present, the limit of ~ 1.56 msec is not very re- 
strictive, and is consistent with basically all realistic models of high density matter. Discovery of a 
sub-millisecond (P ^ 0.5 msec) pulsar would allow to distinguish more directly between competing 
supernuclear equations of state p], 11^, 47|. 



D. Pulsar Glitches 



Pulsar periods are observed to increase gradually with time, implying that rotational energy of the 
star is being lost (primarily by electromagnetic dipole radiation). However, several pulsars have 
been observed to undergo sudden decreases in the rotation period, which are believed to originate 
from the transfer of angular momentum between different parts of the star. The current leading 
model 1 48, 4£] suggests that the source of excess angular momentum is the inner curst, and thus 
the magnitude of glitch phenomena can be used to set a lower limit on the fraction of the total 
moment of inertia the inner crust must carry. Recent results |5^] suggest the crust of the Vela 
pulsar must carry at least 1.4% of the total moment of inertia of the star, which is sufficient to 
rule out equations of state which are very soft in the range 1 < p/pnuc ^ 2. Possible progress 
in theoretical modeling of glitch phenomena could provide further limits on the properties of high 
density matter near the nuclear saturation density. 
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VI. Thermal Properties of Matter in White Dwarfs and Neutron 
Stars 



Fermions in both white dwarfs and neutron stars are expected to have degenerate energies which 
are much larger than their thermal energies. However, compact objects do have finite temperatures 
that are a relic of their progenitors. Since degenerate matter has very long scattering lengths for 
individual particles, the cores of both white dwarfs and neutron stars are expected to be excellent 
thermal conductors, and practically isothermal. 

Finite temperature effects on the properties of condensed matter in compact objects has become 
a very active theoretical and observational field in recent years, especially in the case of white dwarfs. 
A finite temperature leads to cooling through thermal emission processes which are potentially 
observable. Examining the cooling history (i.e., temperature vs. age) of compact objects thus 
provides empirical evidence regarding the thermal properties of high density matter. Here we can 
only briefly mention some of the topics of current interest. 



A. Thermal Effects in White Dwarf Structure 



The region most affected by a finite temperature is the low density outer-layers of a white dwarf, 
which may include a thin atmosphere composed of helium and possibly also hydrogen. For lower 
densities the electrons are not degenerate and a surface temperature of 10^ °K will severely alter the 
equation of state, especially in the regime where electrostatic corrections are important. Modern 



equations of state for a hydrogen/helium atmosphere |51] have been incorporated in several studies 
of white dwarf structure (see ||52| for helium stars, for carbon-oxygen, and fs^ for a recent 
compilation of various compositions). The thermal pressure in the outer layers is generally found 
to be very effective in inflating the stellar radius. As seen from Figure |5|, a surface temperature of 
Tgff = 10^ °K increases the radius by a signiflcant fraction |35, 53] with respect to the cold [ p^ 
models. 

There has been a recent revival of studies of white dwarf thermal evolution, and its dependence 
on the properties of matter at the relevant high densities. We briefly review here the two main 
aspects of these studies, namely cooling and pulsations. 



1. White Dwarf Cooling 

The basic theory of white dwarf cooling was established by Mestel in 1952 |54|. The heat is originally 
stored in the nondegenerate ion-lattice, and the dominant cooling mechanism is photon diffusion to 
the surface from the isothermal core through the nondegenerate outer layers. The simplest results 
are based on Kramer's approximation for the photon opacities in the nondegenerate regime: 

K = , (26) 
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where k is the mean opacity (in cm^gm^^ ) and kq = 4.34x 10^'*Z(l+Xj|/) cm^gm^^ is a composition 
dependent constant, Xh and Z being the mass-fraction of hydrogen and heavy (not hydrogen or 
hehum) elements, respectively. The luminosity through the outer layers, L{r), is dependent on the 
temperature gradient according to the radiative diffusion equation 



(where a = 7.7565 x 10^^^ erg deg~^ is the radiation constant). Combining the above thermal 
temperature gradient with the pressure gradient required for hydrostatic stability (Eq. (1£), see 
Q), one finds that the surface luminosity satisfies 

L = 5,7xl0Vn^^^(|-)T-. (28) 

where fi is the mean molecular weight in the nondegenerate atmosphere and Tc is the core tem- 
perature in the star. For typical white dwarfs (M ~ 1 Mq) with no hydrogen in the atmosphere, 
X = 0, Z K, 0.1, [1 = 1.4, Ye = 0.5, an observed luminosity in the range 10~^ — 10~^ Lq corresponds 
to a central temperature of Tc ~ 10^ — 10^ °K . 

The key elements of Eq. ( [2^ ) are that the surface luminosity (which is observable) is clearly 
related to the mass of the star, the composition and the internal temperature. If the mass is also 
determined separately (as discussed above), luminosity can then provide a direct test of the star's 
thermal and composition profiles. Furthermore, the age of the white dwarf can also be estimated 
based on its cooling time scale, which is dependent on the luminosity function (Eq. the stellar 
mass and the specific heat capacity of a white dwarf material. The latter (in ergs deg~^ gm~^) if 
roughly that of the ion lattice, 

c,(ions) = 2^^ , (29) 
2 fx 

where the factor of 2 comes from the three degrees of collective vibration and rotation, in addition 
to the three single ion degrees of motion. Due to their degeneracy, the specific heat capacity of the 
electrons is suppressed by a factor of ~ ksT/Ep [^], which for a density of 10^ gm cm~^ and an 
internal temperature of ~ 10^ °K is 10"^. 

While models which include more detailed microphysics do reproduce the general results of the 
Mestel model (see, for example, in Refs. |l|, |3^, ^), accurate work that is compatible with 
the quality of current observations must include several important perturbations. Most notably, 
accurate low temperature opacities, semidegenerate electron thermal conductivity and pressure- 
induced ionization must be accounted for in modeling the white dwarf atmosphere. Crystallization 
of the ion gas may also have a significant effect on the thermal history of the star, both as a transient 
source of (latent) heat and by effectively reducing the heat capacity. In particular, below the Debye 
temperature, G^) ~ 4 x lO^p^/^ °K , the specific heat capacity becomes temperature dependent 

= — fe(e^) ■ (30) 

The specific effects of the finite temperature on the thermodynamic properties of matter in white 
dwarf atmosphere are most notable through studies of pulsations (see below). 
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It is noteworthy that for several years there seemed to be a paucity of low-luminosity (and 
therefore, old) observed white dwarfs, which posed a nontrivial puzzle in terms of models for the 
stellar evolution. Only rather recently it was realized that the properties of the finite temper- 
ature atmospheres would make old white dwarfs dimmer than previously anticipated (especially if 
their atmosphere included only a helium layer and no hydrogen). 



2. Pulsation of White Dwarfs 



The basic theory of stellar pulsation is reviewed in detail by Sarbani Basu in this volume. We 
note that studies of pulsations in general allow one to estimate the sound speed profile of the 
star, its temperature gradient, and even the mean molecular weight in its core. It has been well 
established that stellar pulsation provides an effective probe for the properties of matter in 
a star. Commonly referred to as "asteroseismology" , this field is rapidly emerging due significant 
advances in observational instrumentation and theoretical modeling. 

White dwarfs, like many other stars, may undergo stable pulsations which leave a detectable 
imprint on the time dependence of the star's luminosity pO| , |6l| ]. Since pulsations tend to damp 
over time (or the star evolves out of the instability strip, if it is overstable to the pulsations), 
they are most easily observable in younger, and therefore hotter, white dwarfs. Modeling observed 
pulsations of young white dwarfs is particularly important for probing the EoS of semidegenerate 
matter With the aid of theoretical models, studies of white dwarfs pulsations have been 

used for independently restricting the mass-radius relationships and dimensions of the hydrogen 
envelope [61, 64 1, and in general seems to hold much potential for future studies of the details of 



white dwarf physics. For example, some pulsation modes are theoretically predicted to be more 
sensitive to the extent of crystallization of the atmosphere material. Recent analysis of the periods 



and relative magnitudes various pulsation modes in the white dwarf BPM 37093 [63| helps determine 
the extent of crystallization in the atmosphere (about 50% in mass), which compares favorably with 
theory. In another recent example |6^], it is possible to identify, for the first time, the contraction 
rate of the pre- white dwarf star PG 1159-035, by measuring the secular changes over time of the 
periods of various modes. 



B. Thermal Effects in Neutron Star Structure 



Although finite temperature effects are practically negligible when considering the overall structure 
of evolved neutron stars, they are important in assessing their thermal history and as probes of 
matter in the interior of these objects. Some thermal effects on structure do exist in newly-born 
neutron stars (often called "proto- neutron star"). 



1. Proto-Neutron Stars 



Neutron stars are born in supernovae with initial internal temperatures of several times 10 °K and 
initial entropies of ~ 2kB per nucleon (ks is the Boltzmann constant). These values are large enough 
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to impose nontrivial effects on the composition and structure of the neutron star compared to the 
r = case [^]. At such high temperatures the matter is not transparent to the thermal neutrinos 
||l|, and the proto-neutron star must cool through neutrino diffusion to the surface, which 
occurs over a time scale of several seconds (much longer than a dynamical time of milliseconds) 
1 67]. This fundamental prediction was dramatically confirmed by the observed neutrino pulse from 
the supernova 1987A (the closest supernova observed from the Earth in 400 years), which lasted 
for ~ 12 seconds, (see for a review). Further studies of the properties of high temperature, 
super nuclear-density matter are required to model this brief early cooling epoch and its possible 
observational signatures [B^]. 



2. Cooling of Neutron Stars 



After the initial rapid cooling stage which lasts several days, a neutron star settles on a slower 
cooling curve Once neutrinos escape the system freely, they must be continuously produced 

in order to fuel the emission process. Subsequent neutrinos are mostly produced by the so-called 
URCA (/3 and inverse /3 decays as in Eq. (0)) and other processes. The isothermal core at this 
stage is expected to have temperature of several times 10* °K , while at the surface the temperature 
is down to several times 10^ °K - which provides for continuous surface thermal luminosity in soft 
X-rays. This soft X-ray luminosity is much more difficult to detect than the optical ultraviolet 
luminosity white dwarfs. Current observations are generally limited to young (age ^ 10^ yrs) 
neutron stars, where the main cooling mechanism is still the neutrino emission from the core. An 
important consequence of this situation is that estimates of thermal emission from neutron stars 



serve as a probe of the properties of matter in the core |71]. 



Over 20 compact, isolated soft X-ray sources observed with X-ray imaging telescopes have 



been identified as neutron stars |72]. In most cases the source is either too faint to allow for 
a spectroscopic determination of the surface temperature, or the emission is dominated by other 
phenomena (most likely magnetospheric emission). However, in a hand- full of cases surface thermal 
emission has been strong enough to be detected and the surface temperature has been deduced to 



within a factor of two [ 72 1 . The measurements are still not accurate enough to determine whether 
the observed layer is the actual surface or a possible hydrogen atmosphere, which would have 
different emission properties, but they are sufficient to place some constraints on the rates of 
neutrino emitting processes in the core [^ ]. 

Perhaps the most notable conclusion to date of observed neutron star cooling rates is that 
they support the theoretical prediction that the nucleons in the core couple to a superfluid state. 
Theoretical models suggest that the strong interactions will pair the neutrons in the core in a 
^Pi superfluid and the protons to a ^5*0 superconductor with critical temperatures of the order of 
10^ °K |7^. Neutrino emitting processes, such as /3-decays must break a coupled Cooper pair before 
its constituents can participate in the decay. The main effect of nucleon superfluidity is therefore 
damping the efficiency of neutrino emission from the core (there is also a secondary effect due to 
the modulation of the heat capacity, which first increases discontinuously as the star cools to the 
critical temperature, and then decreases exponentially at lower temperatures). Observed neutron 
star surface temperatures are apparently too high to be consistent with the cooling rates predicted 
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for normal core, but they agree with the suppressed coohng rates when the nucleons couple to a 
superfluid state |75, |7^. Further progress in analyses of the matter in the interiors of neutron stars 
is expected through measurements by NASA's recently launched Chandra X-ray satellite. 



VII. Concluding Remarks 

The theoretical and observational study of compact objects remains one of the most exciting fields 
in modern astronomy. In essence, this research is also an exploration of the properties of condensed 
matter at extreme densities. Predictions regarding the properties of white dwarfs and neutron 
stars serve to test our understanding of matter at these high densities, while theories of high 
density matter serve as a basis for interpreting observational results regarding these objects. Most 
exciting, these objects bring together all four of the fundamental forces of nature and probe regimes 
not accessible in the terrestrial laboratory. They provide the most numerous and accessible sample 
of objects where relativistic gravitation - general relativity - plays a role in determining their 
physical properties. 

In this review we have described the tight interconnection between the microscopic (local) 
properties of condensed matter at high densities and the macroscopic (global) properties of white 
dwarfs and neutron stars. While the fundamental principles of cold, high density matter are believed 
to be well understood, and are generally consistent with observations, key questions still remain, 
and new observations may give rise to new puzzles. The current boom in capabilities of Earth- 
bound telescopes and satellite instrumentation promises that many more puzzles - and hopefully, 
answers - are in store regarding the nature of cosmic matter at high densities. 
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